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بهبود بخشی نمودار پارتو در بهینه‌سازی توپولوژی چند هدفه با المان های چندضلعی غیریکنواخت * 


۵ (۱) ِ ۹ )۲( 
مهناز اکبرزاغی بهروز احمدی ندوشن 


چکیده یکی از روش‌های پیشنهادی برای کاهش وزن سازه‌هاء معرفی فضاهای نعالی در سازه است. این موضوع ایده‌ی اولیه شک لگیری 
بویا زین خر ری است, یکی رآ مختکا یقن هش تنل قاری تزر ترش تا ن‌ هی جهازخل ی مهتم سس رخ سخلی آیجاد 
پدیده‌ی شطرنجی است. به‌طو رکل ی گسسته‌سازی که تخحمین بهتری از دامنه پیوسته ارائه دهد ایجاد پدیده‌ی شطرنجی را کاهش می‌دهد. ینک 
زاهفکا یرانق جل وگیر ییاز پدیلوی شطرتج ید مسا بهیه‌سار ی تویرلرژی أسفاهه از لجان‌های جتصلعی است, اف رآین قاله: در بقال بت 
دامنه طراحی محدب و غیر محدب برای بررس ی کاربرد المان‌های چندضلعی غیریکنواحت در بهینه‌سازی توپولوژی چندهدفه تحلیل شده و 


نتایج نشان دهنده کاهش زمان محاسبات و بهبود نمودار پارتو با المان غیریکنوانحت می باشند. 


واژه‌های کلیدی بهینه‌سازی توپولوژی چندهدفه. نمودار پارتو المان‌های چندضلعی. الگوریتم لوید. نمودار ورونوی. 
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* تاریخ دریافت مقاله ۹0/4/۲0 و تاریخ پذیرش آن ٩/۱/۳۰‏ می‌باشد. 


(۱) نویسنده مسئول. کارشناسی ارشد. گروه عمران. دانشکده فنی و مهندسی و دانشگاه یزد 
جصمم. تمصع 6 نم 2 2۲هع1ج. مصطع]۱۷ :1تعصر۴ 


(۲) دانشیار. گروه عمران دانشکده فنی و مهندسی و دانشگاه یزد 
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مقد مه 
درگذشته از المان‌های چهارضلعی و مثلئی درزمینه‌ی 
بهینه‌سازی توپولوژی استفاده می‌کردند. بااین‌حال این 
المان‌ها به طور متداول. دچار ناپایداری عددی ازجمله 
پدیده شطرنجی در مسائل بهینه‌سازی توپولوژی 
یشان که قی ترا از یلرتنک رای اوسن .بردن 
ناپایداری عددی از جمله پدیده شطرنجی استفاده کرد. 
استفاده از فیلترینگ می‌تواند مشکل بهینه‌سازی را 
اقتاتقن دهد [ ی ]زک راهعحل .من تلع از نیع تروق 
بشین مک کاه تم اوه از ام تون که شاف 
(لاهشییی) اسر افظافه او البان‌های. کت ای 
بدون اعمال هرگونه محدودیت. مشکل پدیده شطرنجی 
را از بين می‌برد و علت آن این است که المان‌های 
شش‌ضلعی یا یه یکدیگر اتصال ندارئد یا توسط دو گره 
و یک ضلع با یکدیگر اتصال دارند [3]. یکی از 
مشکلات. استفاده. از شبکه‌های. بکبراشت. مانند 
المان‌های لانه‌زنبوری. مشکل گسسته‌سازی دامنه 
طراحی و دقت در نشان‌دادن بارگذاری و شرایط مرزی 
است. المان‌های چندضلعی غیریکنواخت با فراهم کردن 
گسسته‌سازی انعطاف‌پذیرتر در دامنه‌های پیچیده. 
می‌تواند در گسسته‌سازی بهینه‌سازی توپولوژی مفیدتر 
باشد [4]. 
در ادامه شال تبر 1488 با دامشه و شرایط 
بارگذاری مطابق شکل (۱) بهینه‌سازی توپولوژی را 
بدون محدودیت اضافی با استفاده از روش‌های مختلف 


مش‌بندی بررسی شده ایستا: 


شکل ۱ دامنه و شرایط با گذاری تیر ۷188 با مشخصات مصالح 


استفاده شده در آن 


در شکل ۲( نتایج بهینه‌سازی تویولوژی ارائه شده 


نشریة مهندسی عمران فردوسی 


بهبود بحشی نمودار پارتو در بهینه‌سازی توپولوژی چند هدفه.. 


که ردیف بالا مش‌بندی با المان چهارضلعی و ردیف 
وسط مش‌بندی با المان شش‌ضلعی (لانه‌زنبوری) 
استفاده شده‌است. همانطور که در شکل مشاهده می- 
گردد مش‌بندی با المان‌های چهارضلعی. الگوهای 
شطرنجی ایجاد می‌شود . الگوهای شطرنجی دارای 
سختی کاذب هستند که در مقابل. راه‌حل مش‌بندی‌های 
با المان‌های چندضلعی (لانه‌زنبوری و غیریکنواخت) از 
چنین ناهنجاری‌های آزاد هستند و سختی کاذب ندارند 
[3]. باید به این نکته اشاره کرد که به‌منظور مدل کردن 
یک دامنه با مرز مستقیم استفاده از المان شش‌ضلعی: 
باید یکلایه المان مثلثی و چهارضلعی در امتداد مرز 
قرار دهیم که این در مش‌بندی یکنواخت مناسب نیست 
و از معایب مش‌بندی شش‌ضلعی می‌باشد [5]. 


شکل ۲ راه‌حل توپولوژی برای مسئله تیر 118 بدون هرگونه 
لانه‌زنبوری. پایین: مش بندی جندضلعی غیریکنواخت 


شکل ۳ طرح نهایی تیر 1488 با توجه به طرح میشل [6] 


همان طور که مشاهده شد نتایج بهینه‌سازی 


سال سی و یکم. شمارهُ دو ۱۳۹۷ 


مهنا ز اکبر زاغی- بهروز احمدی ندوشن 


قوریز اوق میس تلم بسا اسان بان سنا خسن 
غیریکنواخت ازالمان‌های مربصی و لانه زنسوری بهشر 
عمل می‌کند و به طرح نهایی ارائه‌شده توسط میشل 
مطابق شکل () نزدیک تر است. بنابراین در ادامه اینن 
مقاله از مش‌بندی غیریکنواخت استفاده شده و مزایای 
آن در مسائل بهینه‌سازی توپولوژی چندهدفه با حل دو 
مثال با دامنه‌ی محدب و غیرمحدب بررسی می‌گردد. 


ساخت المان‌های جندضلعی غیریکنواخت 

المان‌های چندضلعی غیریکنواخت با فراهم کردن 
گسسته‌سازی انعطاف پذیرتر در دامنه‌های پیچیده. 
می‌تواند در گسسته‌سازی بهینه‌سازی توپولوژی مفیدتر 
باشد. در ساخت این المان‌ها. از نمودار ورونوی 
به عنوان وسیله برای تولید المان‌های چندضلعی نامنظم 
ستفاده می‌شود [7]. از ویژگی‌های جالب این روش این 
ست که سطوح تصادفی و با همسانی هندسی از قرار 
دادن تقاط دلخواه و کاملا تصادفی به دست ین ابا و در 
دامه از الگوریتم لوید برای یکنواخت کردن المان‌ها 
ستفاده می‌شود [8], 


مش‌بندی اولیه و دباگرام ورونوی. برای تولید 
مش‌بندی‌ها غیریکنواخت. از نمودار ورونوی برای 
گسسته‌سازی دامنه استفاده می‌شود. گسسته‌سازی دامنه 
با استفاده از روش المان‌های غیریکنواخت دارای تقریب 
مناسبی از مرز است. هم چنین از روش وید برای ایجاد 
یک توزیع دانه‌های یکنواخت و درنتیجه ساخت مش- 
بندی با کیفیت بالا استفاده می‌شود. مش‌بندی 
تا کی رن انم ای مرخ قاط تضادق چی دام 
۵ و همچنین نقاط کمکی برای تقریب شرایط مرزی 
04 ایجاد می‌شود. به‌طورکلی از روش زير برای مش- 
بندی اولیه استفاده می‌شود: 
۱. در داخل دامنه 4 مجموعه نقاط تصادفی به تعداد 
موردنظر تولید می‌شوند. به اين مجموعه نقاط ,۳ 


2 


می‌گویند. 


تال سی و یکم. شمارة دو» ۱۳۹۷ 


۸۹ 


۷ برا ایجاد ظریبت اسب از داست تقاط داغلی بان 3 
نسبت به لبه دامنه منعکس شوند. مجموعه نقاط 
به‌دست آمده را نقاط کمکی میم" می‌نامند. 

۳ دیاگرام ورونوی با نقاط ,۲ لایر < ۳ ساخته 
ی تتر3: 

گسسته سازی چندضلعی از دامنه توسط المان‌های در 
ارتباط با نقاط ایجاد می‌شود. 

این رون در کل( برای تلف ک نان 
پنج ضلعی. استفاده‌شده است [4,9], 


شکل ۶ تولید المان غیریکنواخت. راست: قرارگیری نقاط تصادفی. 
چپ: ساخت دیاگرام ورونوی با مجموعه نقاط تصادفی و کمکی 


(مربع). 


برای ساخت چنین مش‌بندی‌هایی بعد از پراکنده 
درون دامنه, هیچ تأثیری در تقریب دامنه ندارد. بنابراین 
شرط زیر برای انعکاس نقاطی که در نزدیکی مرز لحاظ 
می‌دهد: 
٩( (۱)‏ ,)0 > |(7) رم| 
که در آن ((۳,6) یک مقدار فاصله متناسب با 
عرض المان است که طبق رابطه زیر محاسبه می‌شود: 


2 2 


)۲( ك( 6 <: (() ,)0 


برای پیدا کردن بازتابی از * نسبت به نزدیک‌ترین 
نقطه مرزی استفاده می‌شود. بازتاب‌ها با علامت 


(۳)0 مشخص می‌شود و طبق رابطه زیر به‌دست 


نشريهة مهندسی عمران فردوسی 


۳( 20200۲۷000 << )م1 


که در آن ۷4 جهت نزدیک‌ترین نقطه مرزی را 
می‌دهد و تابع فاصله 800 به‌صورت زیر تعریف 
گردد. اگر 44 زیرمجموعه‌ای از 18 با شرایط مرزی 
مشخص باشد. تابع فاصله مربوط به 4 با ۲ + 4:1 
با رابطة زیر تعریف می‌شود: 


2 7۷۱ 5 ۹ ومء 1۳۳ 50 < )وه 


در آن 040 نشان‌دهن له مسرز 420 اسست؛ 
بنابراین |7۱ - | فاصله اقلیدسی بین دونقطه‌ی ۶ و۷ 
روی مرز دامنه است و تابع م5 در رابطه زیر تعریف‌شده 
است : 


61 
2 


(0) )21-200 60 وه 


در بازتاب نقاط گاهی همانند شکل (۵) نقطه‌ای 
مانند ‏ دارای فاصله‌ی مساوی با دو مرز می‌باشد. آنگاه 
)۷0۵ به‌صورت واحد وجود ندارد و در این حالت * 
باید نسبت به هر دو مرز منعکس شود. 


, 1 
291 

ی( 

ی 


>ذر_" ۱ 7 سس 
٩۱۸ ۰‏ (ام 
ل* «. 


شکل ۵ نقطه در بیش از یک نقطه. فاصله مساوی از 060 است 


همچنین در بازتاب نقاط نسبت به مرز باید توجه 
کرد که انعکاس یک نقطه نسبت به نزدیک‌ترین مرز 
مطایق شکل (1) ممکن است برای یافتن گوشة مناسب 
کافی نباشد و این مشکل توسط انعکاس نقطه نسبت به 


هر دو بخش مرز مجاور, حل‌وفصل می گردد. 


نشریة مهندسی عمران فردوسی 


بهبود بحشی نمودار پارنو در بهینه‌سازی توپولوژی چند هدفه... 


بعد از انعکاس نقاط و با استفاده از نقاط تصادفی, 
تقسیم‌بندی وروئوی صورت گرفته و با استفاده از آن‌ها 


مش‌بندی ایجاد می‌شود. 
():1۲ 
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شکل 1 برای ایجاد یک‌گوشه دقیق نیاز به بازتاب نقاط نسبت به 


هر دو بخثر مرز در گوشه است. 


الگوریتم لوید. نظم نمودار ورونوی به‌ طور کامل 
توسط مجموعه نقاط تولیدشده تعیین می‌شود. یک 
مجموعه مولد تصادفی يا شبه تصادفی ممکن است به 
که کسهستادش کنه ی ام اسفادهی لیا الناق 
محدود مناسب نیست. منجر شود. بنابراین باید یک 
مجموعه‌ای از تقسیم‌بندی ورونوی که یک سطح 
بالاتری از نم را دارد. را بدست آورد. یکی از 
الگوریتم‌های پرطرفدار برای ساخت مش‌بندی 
چندضلعی: طرح تکراری الگوریتم لوید است [8]. 
الگوریتم لوید با کاهش انرژی به صورت محلی همگرا 
می‌شود. مراحل انجام الگوریتم لوید به شرح زیر قابل 
بیان است: 

۱. ساختن نمودار ورونوی در ارتباط با نقاط 

۲. محاسبه مرکز جرم هر سلول 


سال سی و یکم. شمارهُ دو ۱۳۹۷ 


مهناز اکبرزاغی- بهروز احمدی ندوشن 


رفتن به گام بعد مگر اينکه به همگرایی رسیده باشد. 
در الگوریتم لوید. مرکز ثقل هر قسمت به‌صورت 
زیر محاسبه می‌شود و جایگزین نقاط اولیه می‌گردد: 
69۹ ۲ 
۱( سس تاد یوم و < 7 
( ۵ وبا 
در آن (۳ تابع چگالی تعریف‌شده روی 2 است. 
ازاین‌رو در تقسیم‌بندی‌های ایجادشده. هر نقطه‌ی ۷ 
تولیدشده هم‌زمان با نقطه ثقلی ۷۰ نیز مرتبط است. بعد 
از جایگزین شدن نقاط .۷ به‌جای نقاط قبلی. انحراف 
نقاط تولیدشده در هر ناحیه ورونوی. توسط تابع انرژی 
طبق رابطه زیر اندازه‌گیری می‌شود: 
)۷( ۷ 169 ۳ و2 ی )10 ,)۶ 


دیده می‌شود که تابع انرژی در تکرار متوالی از 
الگوریتم وید کاهش می‌یابد که در رابطه زیر نشان 
داده‌شده است : 
)۸( (۴,۸۵)ع > (۵بیبز؟)ء 

به این معنی که الگوریتم لوید می‌تواند به عنوان 
یک روش نزولی برای تابع انرژی مشاهده شود [10]. 
ماه شود که ی بصی با اقرانش کزان الکتروین 
لوید بهبود می‌یابد و المان های یکنواخت‌تری ایجاد 
می‌گردد و این ویژگی در شکل ( تا )٩(‏ به‌حوبی 
نشان داده شده است. همچنین مقدار کاهش تابع انرژی 
تاش اد اتجرافت تقاط ثر تتفیو کی تک اراس مسا 
ذر شکل (۱۰) نشان داده‌شده است و مشاهده می‌شسود 
که مقدار کاهش تابع انرژی در تکرارهای اول بیشتر از 
تکرارهای بعدی آن است [۱۱]. 


4 ۹ 0 
تک و 
اک ارت 


شکل ۷ تکرار اول الگوریتم لوید 20.4192- 5 


سال سی و یکم. شمارة دو. ۱۳۹۷ 


۹۱ 


.۰ دواد 


شکل ۱۰ نمودار کاهش انرژی با توجه به تکرار الگوریتم لوید 


بررسی توابع شکل و تقریب سازی به‌وسيلة 
المان‌های ایزوپارامتریک 
همان‌طور که از نتایج الگوریتم لوید مشخص است 
المان‌های ایجادشده. به‌صورت چندضلعی‌های 
غیریکنواخت هستند و تعداد اضلاع و شکل آن‌ها باهم 
متفاوت است. در اين مقاله برای تقریب‌سازی توابع 
شکل از المان‌های ایزوپارامتریک استفاده می‌شود. منظور 
از الفان انیوبا اهر یکت, این اس که همه محانسانت بر 
مبنای توابع شکل یکسانی انجام می‌گیرد. برای ایین 
منظور توابع شکل بر روی یک المان استاندارد با 
مختصات ٩9‏ (رو,ی6) < 5 تعریف می‌شود که این 
نوع مختصات نوعی مختصات گرانیگاهی است که مرکز 
دستگاه» مرکز جرم آن المان است. فرمول المان محدود 
۷ ضلعی که شامل درون‌یابی توابع شکل و تقریب- 


نشریهة مهندسی عمران فردوسی 
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سازی با المان ایزوپارامتریک بررسی می‌گردد. 


زیر تعریف می‌شوند: 


(8) 


شکل ۱۱ تصویری از یک المان و مثلث‌های ساخته‌شده با آن 


درون‌پابی تابع شکل مربوط به گره ام به‌صورت 
زیر تعریف می‌شود (شکل ۱۱): 


(6) :0 
۱۰ - ۰ 
ٍِ" هر ۳ > (6):ا 
1 
۱۱2( < (0)6 


ریز طرو۳)ه6(۵,رطری- )۸ 


که ۵ سیاحخت:مقلت. با مرقیت زنوس داخل 
زان میدافل 

بدیهی است که توابع شکل در امتداد لبه‌ی 
چندضلعی خطی هستند و خاصیت دلتای کرنوکر را 
برآورده می‌کند و همچنین توابع شکل غیر منفی هستند 
و جمع آن‌ها برابر با یک می‌باشد. 

با توجه به این خواص, توابع شکل می‌تواند برای 
ساخت نقشه‌هىای ایزوپارامتریک برای هر 9 ضلعی 
محدب استفاده شود برای توضیح کامل‌تر این مبحثه 
خوانندگان به مرجع [12] ارجاع می‌گردند. 

استفاده از مش‌بندی چندضلعی در بهینه‌سازی 
توپولوژی چندهدفه 

هدف اصلی در بهینه‌سازی توپولوژی پیدا کردن چیدمان 


بهینه یک سازه در ناحیه معین برای انتقال بارهای وارده 


نشریة مهندسی عمران فردوسی 


بهبود بحشی نمودار پارنو در بهینه‌سازی توپولوژی چند هدفه... 


به نواحی مرزی است. تنها کمیت‌های معین در مسئله‌ی 
بهینه سازی توپولوژی. بارهای وارده. شرایط تکیه گاهی؛ 
حجم سازه و ممکن است محدودیت‌هایی نظیر موقعیت 
و ابعاد نواحی توخالی و توپر نیز مشخص باشد اما 
اندازه» شکل و ارتباط بین اجزای تشکیل‌دهنده‌ی سازه 
مجهول است. بهینه‌سازی توپولوژی درواقع انتخاب 
هم‌زمان توپولوژی (نحوه‌ی ارتباط نواحی تشکیل‌دهنده 
سازه)؛ شکل و اندازه اعضاء تشکیل‌دهنده سازه می‌باشد. 
بهینه‌سازی توپولوژی به دلیل اینکه در مقایسه با دیگر 
روش‌های بهینه‌سازی حجم بیشتری از مصالح را کاهش 
می‌دهد از درجه اهمیت بالایی برخوردار است [1]. 

همان‌طور که قبلاً بیان شد در طراحی توپولوژی به 
دنبال پیدا کردن مکان بهینه برای مواد ایزوتروپیک 
داد‌شده هستیم یعنی باید تقاطی را که از مواد 
7 
همچنین برای سهولت ساخت و اجرا؛ نتایج طراحی 
باید به نحوی باشد که تمام فضای طراحی تقریباً از 
خا ریز توا یل فان 2 ونر 
شیوه برای حل این گونه مسائل (اصلاح مسئله)؛ 
جایگزین کردن مقادیر پیوسته به‌جای مقادیر گسسته و 
سپس استفاده از نوعی جریمه برای میل دادن جواب به 
مقادیر گسسته (۱-۰) است [13]. یکی از روش‌هایی که 
کارا عفر ترا انم وه مسائل تاره سل تسشن 
نسبی جریمه‌شده یا مدل 91۷۳ است که به‌صورت زیر 
بیان می‌شود [13]: 


)۱۲( 1 < ط, ۲۵ [ظ(6م(ع - 1) +ع] < (06رظ 


در اینجا 000 چگالی و تابع متغیّر طراحی است. 
دلیل کته ارم انتفاو‌شسته ان است کته سکالن 
خصوصیات مصالح را بین ۰ و 7 درون‌یابی می‌کند و 
۳ ماتریس سختی اولیه در المان 1 ام است و ۳ همان 
جریمه برای میل دادن جواب به مقادیر گسسته می‌باشد. 
حال برای بررسی پارامتر جریمه و نحوة انتخاب آن یک 
تس شیر 195 ما دنه خلت (فتک ۱۲ سود 
پزوسی قرارگرفته آشست: تر ظ۸4 با ایام ۳۸۱ کهتون 


سال سی و یکم. شمارهُ دو ۱۳۹۷ 


مهناز اکبرزاغی- بهروز احمدی دوشن 


اکثر مقالات بهینه‌سازی توپولوژی با بار وارده ۱ نیوتون 
و مدول یانگ ا نیوتون بر مترمربع و نسبت پواسون 
۳ در نظر گرفته می‌شود [1,3,14,15] و برای اینکه 
بتوان با سایر مقالات و روش‌ها مقایسه نمود در این 
مقاله همین مثال مورد بررسی قرارگرفته است. برای 
مش‌بندی از ۳۰۰۰ المان تنش مسطح چندضلعی 
غیریکنواشعت استفاده شد. ثتایج حاصل از دو روش 
تابع جریمه پیوسته و گسسته در جدول (۱) آورده شده 


اسنتت: 


جدول ۱ مقایسه پارامتر جریمه گسسته و پیوسته در دامنه محدب 


تیر 1/388 
تابع هدف جریمه 
۱۸/۷۵ ۱/۵ 
۱۳۰ [۱-۱/۵] 
۱۹۶ ۲ 
۱۹۳۸۸ [۱-۲] 
1۱/۹۹۳ ۲/۵ 
[۱-۲/۵] 
۲۰۷ ۳ 
۱/۹ [۱-۲] 
۱/۹۳۲ ۳/۵ 
1 [۱-۳/۵] 
۲۱/۹۹۸ 1 
۱۲۱۱۹۶ [۱-۶] 


همان‌طور که از جدول (۱) مشاهده می‌شود نتایج 
به‌دستآمده از پارامتر جریمه پیوسته نسبت به پارامتر 
جریمه گسسته اعداد بهتری ارائه می‌دهد و در نتیجه 
ادامه این مقاله از پارامتر جریمه پیوسته استفاده‌شده 
است. 

در مسائل بهینه‌سازی چندهدفه چند هدف به‌طور 
همزمان بهینه می‌شوند و توابع هدف در مسائل 
بهینه‌سازی چندهدفه با مش‌بندی مربعی معمولاً کمینه 
کردن نرمی و کمینه کردن حجم [16], بیشینه کردن 
مقدار ویژه کمینه کردن حداکثر تغییر مکان زیر 
بار گذاری‌ها [17]. حداکثر کردن انعطاف‌پذیری به 


سال سی و یکم. شماره دی ۱۳۹۷ 


٩۳ 


فان ونان خی او و رد سیب س اتب اوقت 
مشتانل سای تفه تسانا مه ام 
چندهدفه. حل واحد کلی ندارند و بلکه در این گونه 
مسائل جواب مسئله به‌صورت دسته جواب بهینه پارتو 
ارائه می‌شود. در مقاله. دو تابع هدف نرمی و حجم 
ب‌صورت زیر در نظر گرفته می‌شود: 
۱۳( < (ا ,م)[ ع ] 
(۱ 7« | < (۷)۵ < و 

که در آن »؟ تابع هدف نرمی است که به چگٌالی 
(متغییر طراحی) و جابجایی وابسته است و و تابع 
هدف حجم می‌باشد که فقط به چگالی وابسته است. 
برای ردیابی بهینه‌سازی توپولوژی دو هدفه با توابع 
مشخص شده از نرم‌افزار متلب استفاده‌شده و نزدیک به 
بر موش تمه تشن انا ک م۱۳ ۷ شلما 
مربوط به مش‌بندی چندضلعی غیریکنواخت و ۱۹۰ 
خط مربوط به تعریف توابع هدف و ردیابی نمودار 
پارتو در توپولوژی چندهدفه و حدود۵۰ خط مربوط به 
تقرفت دانیدهای مایت خی رش اش نت ان 
نوشتن کدهای مربوطه با دادن هر دامنه طراحی 
(محدب و غیر محدب) با هر شرایط مرزی می‌توان 
نمودار پارتو و نتایج آن را به دست آورد. 

به دامنه () از فضای "18 یک دامنه محدب گفته 
می‌شود هرگاه برای هر دو نقطه از دامنه ۶ ولا از 0 و 
هر 1 > ] > 0 رابطة زیر برقرار باشد. 
(۱۵) 6 ۷ + 1 - 1) 


در ادامه مثال تیر 2188 تیر را بر اساس توابع 
نرمی و حجم موردبررسی قرار می‌گیرد. شرایط 
بارگذاری و تکیه گاهی سازه موردبررسی در شکل (۱۲) 
مشخص‌شده است و نتایج حاصل از بهینه‌سازی 


توپولوژی دو هدفه در جدول (۲) ارائه گردیده است. 


نشريهُ مهندسی عمران فردوسی 
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شکل ۱۲ مسئله تير ۰۱۷1۳813 هندسه و شرایط مرزی با دامنه 
محدب 
همچنین نمودار پارتو مربوط به تیر ۷1813 با شرایط 
بارگذاری توضیح داده شده در شکل (۱۳ ردیابی م ی 
گردد. همان‌طور که مشاهده می‌شود در مسائل با داشتن 


دو تابع هدف. با یک نمودار پارتر و مجموعه‌ای از 


بهبود بخشی نمودار پارنو در بهینه‌سازی توپولوژی چند هدفه... 


توپولوژی‌های مختلف روبرو هستیم. در کسرهای 
حجمی بزرگ. به دلیل مبهم بودن نحوه‌ی جیدمان مواد. 
نمی‌توان از روی آن‌هاء توپولوژی نهایی مناسب را 
برحی از اعضا. از توپولوژی سازه حذف می‌شوند. 
همان‌طور که از روی نمودار دیده می‌شود که در 
محدوده‌ی میانی نمودار نیز تراکم نقاط نسبت به هم 
بیشتر می‌شود و ازاین‌رو جواب مناسب بهینه توپولوژی 


در نمودار پارتو در محدوده وسط نمودار قرار می‌گيرد. 


حدول ۲ برخی از نتایج بهینه‌سازی توپولوژی دو هدفه در تیر ۲188 با ۰ المان 


کسر حجمی: ۱/۸ 


نرمی: ۱۳۵/۱ 


کسر حجمی: ۰/۹۵ 
نرمی: ۱۵۳/۹ 


کسر حجمی: ۰*۱۳ 


نرمی: ۱۹۷۱/۷ 


کسر حجمی: ۰,۱۸۷ 


نرمی: ۳5:۸۵ 


کسر حجمی: ۰,۱۰۷ 


نشریةٌ مهندسی عمران فردوسی 


کسر حجمی: ۰/۸۵ 


توتی ۶ ۱۳۱/۲ 


کسر حجمی: ۱/۷ 


نرمی: ۱۸۵/۹ 


کسر حجمی: ۰/6۵ 
نرمی: ۱۳۹/۳ 


کسر حجمی: 2۰۹,* 


نرمی: ۲۳۹/۰ 


کسر حجمی: ۰,۳۳۸ 


نرمی: ۳۰۰/۵ 


کسر حجمی: ۱/۹ 


نرمی: ۱۳۸/۶ 


کسر حجمی: ۰/۷۵ 
نرمی: ۱۳۹/۳۲ 


کسر حجمی: ۱/۹ 
نرمی: ۱۹۱۳/۷ 


کسر حجمی: ۱" 


نرمی: ۳۰۳/۶ 


کسر حجمی: ۰,۵۳۳ 
نرمی: ۱۸۷ 


سال سی و یکم شماره دی ۱۳۹۷ 


مهنا ز اکبرزاغی- بهروز احمدی ندوشن 


»۵ ۵ اا /ا ا "7۳ ۲ ناسآ 
کسر حجمی 
شکل ۱۳ ردیابی بهینه پارتو تیر 183 با المان‌های چندضلعی 
غیریکنواخت 


برای بررسی اعتبارسنجی, نمودار پارتو به‌دست- 
آمده با نمودار پارتو توپولوژی به‌دست‌آمده از روش 
سوروش [1] با مش‌بندی چهارضلعی (مربعی) و با 
توابع هدف مشابه مقایسه می‌گردد. هر دو مثال برای تیر 
با ۱ اسان س آغ من بغلعن سای ور 
تعداد ۱۰*۲۰ المان برای مش بندی مربعی اجراشده 
است و توابع هدف نرمی و حجم به یک گونه 
تعریف‌فنلهاند. برای:مقایسته نمودارهنا هر دو تموداز 
پارتو برحسب نرمی و کسر حجمی در شکل (۱۶) رسم 
شده است. 


۳0 وو۲20-5 
موه ومورود و 


شکل ۱۶ مقایسه بهینه پارتو تبر ۷131 با المان‌های چندضلعی 


غیریکنواخت و المان‌های مربعی در مقاله سوروش [1]. 


همان‌طور که از شکل (۱۶) مشخص است. با 
مقایسه دو منحنی پارتو دیده‌می‌شود که نتایج بهینه‌سازی 
توپولوژی با مش‌بندی چندضلعی از نتایج به‌دست‌آمده 
از روش مش‌بندی مربعی. بهتر می‌باشد و نمودار پارتو 
به‌دست‌آمده از روش مش‌بندی چندضلعی در پایین 


نشریة مهندسی عمران فردوسی 
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نمودار پارتو با مش‌بندی مربعی قرار می‌گیرد. 
به‌عنوان‌مثال در یک حجم مشخص استفاده از مش‌بندی 
چندضلعی غیریکنواخت. مقدار نرمی کمتر و بهتری را 
در اختیار ما قرار می‌دهد. حال با توجه به مثال فوق که 
ردیابی منحنی پارتو در حالت مش‌بندی باالمان‌های 
جندضلعی به جواب‌های قابل‌قبول‌تری می‌دهد. 


من چنتضلفی قر یکتواخت 
المان مربعی _. ۲ مان ,۴ 
الما ِ 
1 .۱ 
ح ۰۰ ۹ 
۳ 9 
کب ۱۸ 
كت 
د 
- فِ 
<< 
هه 
لز ۳۵ ۳ یذ ۰/۶۵ ۴ ۵ / کف وز 
کسر حجمی 


شکل ۱۵ مقایسه بهینه پارتو تیر ۷1198 با المان‌های مربعی و 
المان‌های جندضلعی غیریکنواخت 

در ادامه تأثیر تعداد المان‌ها بر روی نتایج منحنی 
پارتو بررسی می‌گردد. به این منظور مسئله تير ۷۳85 را 
با تعداد المان‌های غیریکنواخت ۸۰۰۰۶۰۰ و ۱۲۰۰ را 
اجرا کرده و نمودار پارتو را برای هر یک از آن‌هابه 
دست آورده و با نتیجه منحنی پارتو با المان مربعی با 
تعداد ۱۰:۲۰ مقایسه می‌شود. 

در هربار اجرای برنامه‌ی متلب. زمان را با دفت 
دقیقه با رایانه شخصی با مشخصات ( 0۳1۱:00۲2 
60 2.206 76670 112) اندازه‌ گیری نموده و 
تأثیر استفاده از المان‌های چندضلعی غیریکنواخت و 
همچنین افزایش تعداد المان بر روی زمان اجرای برنامه 


موردیررسی قرار داده شد. 
جدول ۳ مقایسه زمان اجرای برنامه متلب برای المان‌های 


چندضلعی غیریکنواخت و مربعی 
درصد کاهش ۳ المان تعداد 
زمان نت غیریکنواخت المان 
02 ۰ دقبقه ۱ دقیقه 33 
۱ ۱ دقبقه ۵ دقیقه ۸.۰ 
1۳ ۷ دقیقه ۸ دقیقه ۱۳.۰ 


سال سی و یکم. شمارهُ دو ۱۳۹۷ 


۹۹ 


خسن قرو که از بجتزل شفنب 
افرایش تعداد المان‌ها و ویز شندن مش‌بندی‌هاء زمان 
محاسبات افزایش می‌یابد. همین‌طور که در شکل (۱۵) 
دیده می‌شود ریزترشدن مش‌بندی‌هاء نتایج حاصل از 
منحنی پارتو را بهبود داده و دقت محاسبات بیشتر می- 
شود. زمان اجرای برنامه در مش‌بندی چندضلعی غیر- 
یکنواخت از ثلث زمان اجرا در مش‌بندی مربعی کمتر 
است و درنتیجه می‌توان در زمان کمتر به نتایج 
قابل قبول‌تری با المان‌های غیریکنواخت دست‌یافت. 

در طراحی‌های گسترده و سازه‌هاء افزایش تعداد 
المان در بعضی موارد مقرون‌به‌صرفه نمی‌باشد و باید 
جنبه اقتصادی و زمان اجرا را در نظر گرفت. افزایش 
تعداد المان در چندضلعی غیریکنواخت. در نمودار 
پارتو تغییرات زیادی ایجاد نمی‌کند. درنتیجه با مش- 
بندی چندضلعی غیریکنواخت می‌توان با تعداد المان 
کمت جواب‌های قابل قبولی را به دست آورد. 

طبق نتایج ارائه‌شده مشاهده می‌شود که زمان 
محاسبات. بیش از 1۵ درصد کاهش می‌پابد. علت 
کاهش زمان محاسبات را می‌توان کم شدن زمان 
محاسبات مرتبط با فیلترینگ دانست. درنتیجه چون در 
اکثر طراحی‌ها تابع هدف کمینه کردن وزن یا حجم 
بجای کمینه کردن هزینه استفاده می‌شود با استفاده از 
بهینه‌سازی توپولوژی چندهدفه با المان چندضلعی 
غیریکنواخت. می‌توان سازه‌هایی سبک وزن. کم قیمت 
با بازدهی بالا طراحی نمود. کاهش زیاد زمان محاسبات 
از یک سو و نزدیک بودن طرح نهایی تیر ظظ/۷ 
ارائه‌شده توسط المان‌های غیریکنواخت. به طرح بهینه 
ارائه‌شده توسط میشل [6] از سو دیگر علت‌های اصلی 
اقتصادی بودن روش چندضلعی غیریکنواخت می‌باشد. 

در ادامه» مسئله معروف تير 5010001106 با دامنه 
طراحی و شرایط تکیه‌گاهی مانند شکل (۱7) را با تعداد 
۰ المان به‌عنوان نماینده‌ای از مسائل توپولوژی 
دامنه‌ی غیرمحدب مورد بررسی قرار می گیرد. 
مایت این عال» ماد سکعصانت: اراتهشده خر 
ساير مقالات برای این مسئله. بار وارده ۱ نیوتون و 


مدول یانگ ۱ نیوتون‌برمترمربع و نسبت پواسون ۰/۳ در 
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بهبود بحشی نمودار پارنو در بهینه‌سازی توپولوژی چند هدفه... 


نظر گرفته‌شده است. 

جدول (4) نتایج بهینه‌سازی ارائه‌شده در برنامة 
متلب را نشان می‌دهد و همچنین نمودار پارتو تیر 
6 با المان‌های جندضلعی غیریکنواخت در 
شکل (۱۷) نشان داده‌شده است و همچنین تأثیر افزایش 
تفداد المانزیی وق امیش عحلب درک 9 


موردبررسی قرارگرفته است. 


شکل ۱۱ مسئله تیر 0061106 هندسه و شرایط مرزی با دامنه 


جیز معلی 


۰/۳۵ ۰ ۰/۴۵6 ۰ ۰/۵6 ۴ ۰/6 ۰ ۰/۷5 ۰۸ ۰۸8 ٩ 
کر حجمی‎ 


شکل ۱۷ ردیابی بهینه پارتو تير ۹60601106 با ۳۰۰۰ المان 


چندضلعی غیریکنواخت 


شکل ۱۸ بررسی روند تغییر در منحنی پارتو تير 560601106 با 


افزایش تعداد المان 


سال سی و یکم. شمارهُ دو؛ ۱۳۹۷ 


مهنا ز اکبرزاغی- بهرو ز احمدی دوشن 


جدول ۶ برخی از نتایج بهینه‌سازی تویولوژی دو هدفه تیر ۹67060116 با ۲۰۰۰ المان جندضلعی غیریکنواعت 


کسر حجمی: ۰/۸ 
نرمی : ۳۸/۸/۸ 


کسر حجمی: ۰/۸۵ 
نرمی: ۳/۷/۸۹ 


کسر حجمی: ۱/٩‏ 


نرمی: ۳۹۹/۰ 


کسر حجمی: ۰۱/۰ 


برفی ۰۰ ۲۳۹/۲ 


کسر حجمی: ۰,۱۱۹ 
نرمی: ۶۱۳/۸ 


1 
1 


کسر حجمی: ۰*2۸ 
نرمی : 9۳1/۵ 


1 


کسر حجمی: ۰۱۳۸ 


ترمی ۱۳۵ 


1 
1 


کسر حجمی: ۰۱۹۸ 
نرمی: ۱۳۸۹ 
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کسر حجمی: ۸۱/۰۷/۵" 
نرمی: ۳۱/۶ 


1 


کسر حجمی: ۰,5۳۳ 


نرمی: ۳۷/۰ 


کسر حجمی: ۰:2۸ 


نرمی: 3/۸ 


1 
1 


کسر حجمی: ۰,۳۸ 


نرمی: 11۵/۰ 


کسر حجمی: ۱/۲۱۰۸ 


ثرمی: ۱۰۱۱/۲ 


کسر حجمی: ۱/۷۵ 


نرمی: ۳۰۳/۸ 
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نرمی: ۳۵/۷ 


1 


کسر حجمی: ۱۰۰۳۰,"* 
نرمی: 16۸/۹ 


کسر حجمی: ۱/۶۱۰۸ 


نرمی: ۶/۶ 


1 


کسر حجمی: ۰۱۳۳۸ 


نرمی: ۸3 


۹۷ 


نشريهة مهندسی عمران فردوسی 


۹۸ 


با توجه به نتیجه به دست‌آمده برای تیر 
6 مشاهده می‌شود که می‌توان از مش‌بندی 
غیریکنواخت برای به دست آوردن طرح توپولوژی 
برای دامنه‌های غیرمحدب استفاده کرد و نمودار پارتو را 
ردیابی نمود. با استفاده از المان‌های غیریکنواخت می 
توان به‌راحتی هر دامنه طراحی غیرمحدب و یا دارای 
بازشو را مدل کرد و نتایج را بررسی نمود. یکی از 
مشکلات استفاده از شبکه‌های یکنواخت. مشکل 
گسسته‌سازی دامنه طراحی و دقت در نشان‌دادن 
غیریکنواخت. با فراهم کردن گسسته‌سازی انعطاف- 
پذیرتر در دامنه‌های پیچیده می‌تواند در گسسته‌سازی 
بهینه‌سازی توپولوژی کارآمدتر باشد. همچنین با توجه 
بتک (۱۸) امد می‌شره که اف انش سای انیان 
بر روی نتایج تأثیر زیادی ندارد و می‌توان برای دامنة 
غیرمحدب نیز با تعداد المان کمتر به نتایج قابل‌قبولی 


دست‌یافت. 


نتیجه گیری 
با توجه به اينکه روزبه‌روز بر اهمیت بهینه‌سازی در 


سازه‌ها افزوده می‌شود. مهندسین سازه به دنبال راه- 


مرچ 


بهبود بخشی نمودار پارتو در بهینه‌سازی توپولوژی چند هلفه... 


کارهای مناسب و اقتصادی برای یافتن سازه بهینه 
هستند. استفاده از روش مش‌بندی با المان‌های چند- 
ضلعی غیریکنواخت این امکان را ایجاد می‌کند که هر 
نوع سازه با هر شکل و شرایط تکیه‌گاهی را بتوان 
به‌راحتی مدل نموده و مش بندی و سپس بهینه کرد. با 
استفاده از بررسی‌های مختلفی که صورت گرفته است. 
نتایج بهینه‌سازی توپولوژی با استفاده از این مش‌بندی 
بسیار مناسب‌تر از مش‌بندی‌های مربعی عمل نموده و 
دیگر به فیلترینگ که یک روش غیراقتصادی است؛ 
نیازی نیست. مش‌بندی غیریکنواخت به دلیل کاهش 
زمان اجرای برنامه تا بیش از یک‌سوم ((2 
مقرون‌به‌صرفه است و با مش‌بندی چندضلعی 
غیریکنواخت می‌توان با تعداد المان کمتر. جواب‌های 
قابل قبول‌تری را به دست آورد. درنتیجه استفاده از 
بهینه‌سازی توپولوژی چندهدفه با مش‌بندی 
غیریکنواخت. راه‌حل‌های بهتر و دقیق‌تری نسبت به 
مش‌بندی مربعی در اختیار قرار می‌دهد و نمودار پارتو 
بهبود یافته‌ای را در اختیار کاربران قرار می‌دهد. 


مامتا رجمتاهملصتاون زعماممما طا عصتمهه لقصتامن۳۵۲۵۲۵-۵ 10۲ عم ۱۷۸/120 ص199 رنه یک رطفعتناه 
2010(۰) ,665-679 .00 و5 ٩0۰‏ ,42 ۷۵۱ مطم0)ه2تصتاو۵ فکمه ۱۷۲۵6 240 

جر ماه ه تصمتاعمتصتاوه . بزعمآ۵و۵] طاً وم111)1طامافصد. اهممم اظ ‏ رجمعوهاع۳ .زر فصه .۵ راطع 
امسطمیتد. ‏ مقصصتصتصظ.. لهمما.. مه ومزم‌صم0مممم0 دمص رقلته0صامعاممهان . طد . مصتلوع . م2۳۵۵ 
1998(۰) ,68-75 .00 و1 ,0 .16 ,۷۵۱ وطمتاهمتصتامن 

۲۷۷20۵۵۵۵۵۹۹ مصوت ممتاهمتصتامن ‏ بزوماوو۵] ۳ بت مق 2240 بت ومصتلناهط ون بتطه‌علله 1 
۰ :973 ,۷۵۱ ,۱216011215 22060 باامصمتمص مج ماهموتاب۲ ر"واحمصماه لقطمعم‌هط ۳ ممتاها1۳6020 
2006(۰) ,309-316 00۰ و1 

7 ۲۵ فاممصصماه فاص لقصمعبوامط ر.۲ .1 رو۱۵۵۵26 20 بظ۸ رهت۵۳61ظ وبا .6 ومصئلناهط رن بتطه‌علاه ]1 
۰ 6۱9۱۱66۲۱۱9۱ 1۱ ۱۱۵۱۵05 ۱۵۱۵۲۱۵۵ ۵۲ آمیاوز ۱۱66۲۳۵۵۵0۳8۵1 رصع ته هم عصا نصا ه :صماه2تصتانن 
2010(۰) ,671-698 .00 ,6 0۰ ,82 

امسطمنتطه 10۲ عاحمصماه مات قومتم‌وطمع ۲۷۷ امجممبمعصما۳" وبالت) را 2۳0 ,0.۳ رمصتلنتهط بت بتطه‌علاه ]1 
,569-583 .00 و6 ۱0۰ و37 ۷۵۱ مطملاه2تصاون) تام کل ۱۱ مج منطو .صمتمملصتاوی رعمآممم] 
2009(۰) 
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مهناز اکبرزاغی- بهروز احمدی دوشن 


و "00 اون عماممما اهامای ۵۶ فمطاعصرر 0مصفناصاهاوی ۵۶ زو ماه ۸ ویآ.ت) ,۲۹۵7۷2 
(2009) ,217-237 .00 و3 0۰ و37 ,۷۵۱ مطملاه2تصامن ومصتامامونل تاای۲ ۱ مج فتاه 
مجح عصمتای‌ناووه :عممتاهااهعععا تمصمعم۲۷ همهم را رععتاطام2صتت ‏ 20 ...۷۰ وتع۳20 و.م0 وتا 
٩۱۵۱ ۲۵۷۱۵۷, ۷۵ 41, ٩0. 4, 00. 637-676, )1999(‏ ,"عصطاتتمواه 
عمج مومموام-اهمصمم ه مطوم۱ ۰۳۵۱۲ و.ظ ب1 رک۱۷۵۵۵26 20 بخ رتم۲۵۳ .6 ممصلناهط وت بتطه‌علاه ]1 
م20108صصناو۵) منامام علل ۱ صرح امتمتنه رد طما۱۷]2 صا منز عاصمصصماه لقصمع‌بوامم تم تماهتمطهع 
.(2012) ,309-328 .00 و3 0۰ ,45 ۲۷۵1 
احعصعاه ماتص اقصمعبامم 10۲ قصصحممنه زممم۷۵ مص2لصانن ری رطه‌عافظ 24 .۲ رحعتالظ ولا 81626 
335-350 .00 رته‌عصاتوی رما0)20صنام۲ مصنطومم اقصمتمصعماص 19 معط ۵۶ فعصتلععع۵ظ صا و 6طمااهانام‌طمم 
2010(۰) 
همه مصتانام‌جمن ۶۵۲ مصطاتتمعاه 1۱070 عطا ۵۶ ممممعزه7 ۵9" وبا وتا 240 و رماصه‌صعتامصصط ریم پتار 
.(2006) 102-119 .00 ,1 .۱0 44 ,۷۵۱ بکتمصه ۵۱۱۵۱6۵ ۵ آعیامز مر و" عصمتاهاآهععما تمصمع0 ۲۷ 

۱ اکرژاغیفهتان مپتسازی ترپولوژی خسدمدفه در ساژه‌های پتوسته (بایان‌نامة کارشناسی ازشد) دانشگاه بزه (۱۳۹۶). 
۲ 0۳۳۱۵ ۱۱۵۲۱۵۵۵۴۵۱ و قاصمججماه فاص آمصمعامص مصتطم] من ره رتمهتنها2 ۲ فص .از رتهمصت‌لناه 
(2004) ,2045-2066 .00 ر12 ٩0.‏ و6 ۷۵۱ موه ۱۱ م۷ آمعز ۷۱۵۱۱۱۵۲ 
۲ ۱۵۵0۱0۵۲ ماه که ممتاحملصتامن وتامعزها۱۷۲۵۱-۵ بط رقطاظ 0صج .۲ روعناع۵0۲1 رق.ل م1۷]20611 
(2005) و21-26 .00 و1 .۱۱0 ,36 ۷۵1۰ 50107۷270 عصهمصلعطظ ما دعمصج ۸ رت عصصطاتمعاه متاعصوع 
قح متنمنااگ. رد طماه)۱۷ صا مه قمع ممااعتصتامن . عماجم فصن 99 ۳۸ ریم راطع 
2001(۰) 000۰120-1127 و2 0۰ و21 ,۷۵۱ مطمتاهعتصتامه توتفصتآمژه ول بانط 
۲ ۰۳۶۶۱۵1۵0 رم رصتاطعله 20 ریگ .ظ 1222707 و رقاهمه7عط0 ربق مفه‌عتاقلن ر.ظ ب0۲028501 هر 
مط0ه2لصتاون) ومصتامرآم ول ای قمع تمتتعناو همع ۵ عفمصاً 8 عصلونه ظصا ]۱۷۸ ما ممتاه‌متصتامن 
2011(۰) ,1-16 .00 و1 .0 ,43 ,۷۵01 
موجه و0مطامجظط وممط 1 تصمتامتتصاون مماممن]۳ ریم اوه مه ۲۲۰ ر2050عظ 
2003(۰) مطلاظ رهم0۳1٩‏ 
اوه متامعزمام-ا ند هه ۵۶ ممتاهم‌نامید ما رتعنهمممطمگ مضه وم ومعلصعلتامک ریا ۲۱2۲02 
م6۵ 0002۲6 08 6۵06۲۵۵6 ۱۱۵۲۳۱۵۵۵۵ ۳۱ 1 و صعنوع اون زعمآ۵و1۵ ما مصطاتعمع21 
2002(۰) ,۱۱۵۳۵۵6۵۵ 8 ۲واععل 
عصوا صعلوع حصوتصعدام‌عمطر صهتآم‌جمن" ویک رکاما۸۵061-۸ 200 و.۷ رعصقط2 ورگ وقطه۷ وبا مطعط و2 رما 
مج مممای ‏ رعمطمنته . صتتطتاومع .۵۲ مصفامو .. صمتاهملصتاوه .. زعمامجم . 7۵ا60[ما0-)اناظ 


.(2005) ,142-15 .00 و2 0۰ ,30 ۷۵۱۰ مطمللهمنصاون0 وعصتامم وال ات۱۷۲ 
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سال سی و یکم شمار؛ دو ۱۳۹۷ نشریة مهندسی عمران فردوسی 


۱۰۰ بهبود بعشی نمودار پارتو در بهینه‌سازی توپولوژی چند هدفه ... 


نشریة مهندسی عمران فردوسی سال سی و یکم شمارهٌ دو. ۱۳۹۷ 


